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Monolayer transition-metal dichalcogenides possess a pair of degenerate helical valleys in the band 
structure that exhibit fascinating optical valley polarization. Optical valley polarization, however, is 
limited by carrier lifetimes of these materials. Lifting the valley degeneracy is therefore an attractive 
route for achieving valley polarization. It is very challenging to achieve appreciable valley degeneracy 
splitting with applied magnetic field. We propose a strategy to create giant splitting of the valley 
degeneracy by proximity-induced Zeeman effect. As a demonstration, our first principles calculations 
of monolayer MoTe 2 on a EuO substrate show that valley splitting over 300 meV can be generated. 
The proximity coupling also makes interband transition energies valley dependent, enabling valley 
selection by optical frequency tuning in addition to circular polarization. The valley splitting in the 
heterostructure is also continuously tunable by rotating substrate magnetization. The giant and 
tunable valley splitting adds a readily accessible dimension to the valley-spin physics with rich and 
interesting experimental consequences, and offers a practical avenue for exploring device paradigms 
based on the intrinsic degrees of freedom of electrons. 


Introduction . — As promising valleytronics[l- 
7] materials, monolayer transition-metal dichalco¬ 
genides, MX 2 (M = Mo, W and X=S, Se and Te), 
have attracted significant research interest in re¬ 
cent years [8-11]. The fascinating optoelectronic 
properties arise from a pair of degenerate but in¬ 
equivalent valleys in the vicinities of the vertices of 
the hexagonal Brillouin zone (BZ). The principal 
challenge in exploiting the valley degree of free¬ 
dom in electronic and optoelectronic applications 
is the ability to generate and control the valley po¬ 
larization. Dynamical polarization of monolayer 
MX 2 has been achieved by circularly polarized op¬ 
tical pumping, [8-11] based on valley-selective cir¬ 
cular dichroism; that is, the interband transitions 
at two valleys exclusively absorb the left- and right- 
polarized photon, respectively [4, 8, 9]. The sub¬ 
sequent Berry curvature-induced valley Hall effect 
of M 0 S 2 has also been observed in a recent experi¬ 
ment. [12] These findings offer a new paradigm for 
novel optoelectronic devices. 

Dynamical valley polarization based on opti¬ 
cal pumping relies upon transient non-equilibrium 
photo-carrier distribution, and is ultimately lim¬ 
ited by carrier lifetime. Although the fact that 
photo-carriers are rather short-lived for MX 2 (see 
[13, 14] and references therein) is more an issue 
with material quality than a fundamental limit, it 
does present a challenge to utilizing valley polar¬ 
ization experimentally, especially in beyond-proof- 
of-concept devices. An alternative route to val¬ 
ley polarization is lifting the valley degeneracy by 
breaking time-reversal symmetry. Valley degener¬ 
acy splitting will make possible equilibrium val¬ 
ley polarization via conventional doping, although 
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schemes for creating significant valley splitting has 
yet to be developed. Valley splitting in MX 2 has 
been assessed in a few recent experiments, which 
showed that only small valley splitting, 0.1-0.2 
meV/tesla, can be generated by an external mag¬ 
netic field. [15-18] Here, we propose a strategy to 
create giant valley splitting by proximity-induced 
magnetic interactions. We use first principles cal¬ 
culation to demonstrate that in molybdenum ditel- 
luride (MoTe 2 ) on a europium oxide (EuO) sub¬ 
strate valley splitting larger than 300 meV can be 
generated by an induced Zeeman field. Additional 
Rashba field induced by the substrate and Zeeman 
field further splits the degeneracy in the interband 
optical transition between like spins. Based on the 
lifting of valley and spin degeneracies, we discuss 
new avenues for exploring novel spin and valley 
physics by equilibrium doping or dynamical opti¬ 
cal excitation. We also demonstrate that the val¬ 
ley splitting resulting from the proposed strategy 
is highly tunable. 

MoTe 2 /EuO heterostructure . — The structure of 
monolayer MX 2 is shown in Fig. 1(a), which mim¬ 
ics a two-dimensional honeycomb lattice. The cor¬ 
responding BZ and high symmetry points are de¬ 
picted in Fig. 1(b). The inequivalent K± are 
the vertices of the hexagonal BZ, where the di¬ 
rect band gaps are located. Typical low-energy 
band structure of free-standing monolayer MX 2 
is shown in Fig. 1(c). In the absence of inver¬ 
sion symmetry, the bands are not spin degener¬ 
ate and the electronic states with up/down spins 
are represented by blue/red lines with up arrows 
(t) /down arrows (1). We now introduce a first en¬ 
ergy scale pertinent to the important valley-spin 
physics. The spin splitting within a valley (see 
Fig- 1(c)), arising from 

spin-orbit coupling (SOC). Here, v and c refer re¬ 
spectively to the valence and conduction bands. 
Valleys, K±, are addressed by the index r = ±1. 
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E is the corresponding energy level. When the 
system is invariant under time reversal, the spin 
splitting is opposite for two valleys. Because dif¬ 
ferent d-orbitals dominate the valence and conduc¬ 
tion bands, [ 8 ] the magnitude of the spin splitting 
of the valence bands (~ 150-500 meV) is typically 
larger than that of the conduction bands 3- 
60meV). [19] The relative signs of for the 

conduction and valence bands depends only on the 
type of the metal atom: the signs are opposite for 
M 0 X 2 (as is the case in Fig. 1, where A^p^ is set 
to be positive), while the signs are the same for 
WX2. 



FIG. 1: Atomic and schematic band structures of 
monolayer MX 2 . (a) Top view of a monolayer MX 2 , 
which mimics a bipartite honeycomb lattice. The two 
lattice sites of the honeycomb lattice are occupied by M 
(pink spheres) and a pair of X (yellow spheres). Each 
M is coordinated to six X atoms forming a trigonal 
prism (pink triangles), (b) The Brillouin zone and high 
symmetry points, (c) and (d) Schematic band struc¬ 
ture at the two valleys of monolayer MX 2 without/with 
a Zeeman field, respectively. Blue bands with blue up 
arrows correspond to the up spin states, and red bands 
with red down arrows the down spin states. 

The valley degeneracy can be quantified by an¬ 
other energy scale, 

which is zero in the presence of time-reversal sym¬ 
metry. This important valley degeneracy is funda¬ 
mental to the notion of valley degree of freedom 
and valleytronics based on dynamical valley polar¬ 
ization with chiral optical pumping induced non¬ 
equilibrium carrier distribution [1, 4, 8 , 9]. It is on 
the other hand also highly desirable to break the 
valley degeneracy to exploit the valley and spin 
degrees of freedom, with equilibrium doping or dy¬ 
namical optical excitation, or the combination of 
both. [5] The application of external magnetic field 
breaks time-reversal symmetry and lifts the val¬ 
ley degeneracy, whereby 7 ^ 0 , as exempli¬ 

fied in Fig. 1(d). Systematic measurements re¬ 
vealed that the Zeeman splitting of valley degener¬ 


acy only amounts to ~ 0.1 meV/T. [15] Evidently, 
it is of importance to develop alternative strategies 
to achieve large valley splitting. 

A promising approach is to employ the mag¬ 
netic interaction to lift the valley degeneracy in 
heterostructure composed of monolayer MX 2 and 
an insulating ferromagnetic substrate. In the fol¬ 
lowings, by first-principles density functional the¬ 
ory (DFT) calculations, we take MoTe 2 /EuO(lll) 
as an example to demonstrate a giant valley 
splitting via proximity-induced magnetic interac¬ 
tion. In our first-principles Density functional 
theory calculations, the Perdew-Burke-Ernzerhof 
exchange-correlation functional and the projector- 
augmented wave potentials are used, as imple¬ 
mented in the VASP code. [20] An energy cutoff of 
400 eV for the plane-wave basis is adopted. We use 
al5xl5xl5 Monkhorst-Pack -mesh for bulk EuO, 
while a 15 X 15 X 1 fc-mesh for the pristine MoTe 2 
monolayer and MoTe 2 /EuO. For MoTe 2 monolayer 
and MoTe 2 /EuO, a vacuum slab of 20 A is in¬ 
serted to minimize the interaction between the slab 
model and its periodic images. Atomic coordinates 
are optimized with a convergence threshold of 0.01 
eV/A on the interatomic forces. We use GGA-|-t/ 
calculation to account for on-site Coulomb interac¬ 
tions. [21] The Coulomb and exchange parameters, 
U and J, are set to 8.3 eV and 0.77 eV for Eu 4/ or¬ 
bitals, respectively, while U = 4.6 eV and J = 1.2 
eV for O 2p orbitals. [22, 23] These values give 
a good agreement with the experimental results 
of EuO. [24, 25] The Kohn-Sham Bloch wavefunc- 
tions on the discretized fc-mesh are extracted to 
calculation the optical interband transition matrix 
elements [5, 8 ] and Berry curvature [26]. 

We choose MoTe 2 /EuO heterostructure as an 
example, based on two considerations. Eirst, EuO 
is a ferromagnetic semiconductor with a large band 
gap of more than 1 eV, and offers exchange in¬ 
teraction with ~ 7 /tb spin moment on each Eu 
ions. [27, 28] Second, EuO(lll) single-crystal film 
has very recently been grown on Ir(lll) [29], and 
the lattice mismatch between MoTe 2 and EuO 
(111) substrate is only 2.7%, a reasonable value 
for a commensurate heterostructure composed of 
slightly strained MoTe 2 monolayer on EuO(lll). 
The calculated lattice constant of cubic EuO is 
5.188 A and in-plane lattice constant of EuO (111) 
substrate is therefore set to 3.66 A, which is very 
close to the calculated lattice constant of free¬ 
standing MoTe 2 monolayer, 3.56 A. 

We therefore construct the MoTe 2 /EuO het¬ 
erostructure with a MoTe 2 monolayer placed on 
the Eu-terminated surface of EuO (111) substrate 
composed of 12 Eu/0 atomic layers, as shown in 
Eig. 2(a). The oxygen-terminated surface of EuO 
substrate is saturated by hydrogen, to model a 
semi-infinite single-crystalline EuO or EuO film 
grown on another substrate [29]. Structural relax¬ 
ation reveals a few stable configurations for this 
MoTe 2 /EuO heterostructure, corresponding to rel¬ 
ative shifts of the two materials along the ( 111 ) 
plane of EuO. As detailed in the Supplementary 
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Information (SI), we find that one of these con¬ 
figurations to be substantially more stable than 
the others, which will be focused upon in subse¬ 
quent discussions. In this configuration. Mo atoms 
are located directly on top of Eu ions, favouring 
proximity-induced magnetic effects (Figs. 2 (a) 
and (b)). 



(b) (c) (d) 


FIG. 2 : Atomic and band structures of MoTe2/EuO 
system, (a) Side view and (b) top view of MoTe2 mono- 
layer on EuO ( 111 ) substrate, (c) and (d) Band struc¬ 
tures of MoTe2/EuO with EuO magnetized upward 
and downward, respectively. We use fat band repre¬ 
sentation to indicate the projected weights on MoTe2 
(blue and red) and EuO (gray). The MoTe2 projec¬ 
tions are colored according to the direction of spin mag¬ 
netic moment near K±: blue for up spin and red for 
down spin. The energy scale is zeroed to the Fermi 
level. 

Fig. 2(c) and (d) show the calculated band 
structures of MoTe 2 /EuO heterojunction, with 
EuO magnetized in the upward and downward 
directions, respectively (see insets). The fat- 
band representation corresponds to the projected 
weights of each band onto MoTe 2 (blue and red 
bands), and EuO (gray bands). We see that a few 
bands arise primarily from MoTe 2 in the energy 
range between -1.5 and 0.5 eV, where EuO sub¬ 
strate has only minor contribution. Viewing the 
MoTe 2 bands only, there is a well-defined global 
gap in the energy range between -0.9 and -0.4 eV. 
The direct gaps at K± (~ 0.6 —0.7 eV) indeed cor¬ 
responds to valleys of native free-standing mono- 
layer MoTe 2 . This is also supported by the com¬ 
puted optical selectivity and Berry curvature, to 
be presented shortly. 

The identification of MoTe 2 bands largely free 
of hybridization with the substrate leads immedi¬ 
ately to the key observation that the valley de¬ 
generacy of MoTe 2 bands is substantially lifted. 
The MoTe 2 bands near the gap can be classified 
as spin up and down, respectively colored in blue 
and red in Fig. 2 (c) and (d), as the projected 
spin magnetic moments near Kj- are dominantly 
out-of-plane, along the [Ill] direction of EuO (see 
Fig. S2 of SI). Therefore, the valley splitting can 
be quantified by the magnitude of which 

are in the range from 321 to 419 meV (see Table I). 
Moreover, the smallest energies for the band edge 


TABLE I: Important energy scales of valley and spin, 
for MoTe2 monolayer and MoTe2/EuO heterostruc¬ 
ture. 


Valley splitting (meV) 

a:;! 

a:;| 

Ar 

Ae;7 

MoTez^' (DFT & LEH) 

0 

0 

0 

0 

MoTez/EuO (DFT) 

-342 

-386 

-321 

-419 

MoTez/EuO (LEH) 

-319 

-412 

-340 

-391 

Spin splitting (meV) 

A“’+ 

spin 

A“’+ 

spin 

A"’7 

spin 

A^pL 

MoTez^' (DFT & LEH) 

214 

-27 

-214 

27 

MoTez/EuO (DFT) 

-168 

-449 

-496 

-356 

MoTez/EuO (LEH) 

-142 

-455 

-517 

-348 


1 For free-standing MoTe 2 monolayer, DFT 
calculation and LEH give the same results. 


vertical optical transition without spin flip become 
unequal, with = 886 meV and 

A-pt = E^- - El’- = 930 meV (Fig. 2 (c)), 
and the difference A^jp^ — A^jp^ reaches a substan¬ 
tial value 44 meV, equivalent to the splitting by a 
440-tesla magnetic field. 

Proximity-induced magnetic interactions .— It is 
evident that the substrate induces substantial 
magnetic interaction by proximity effect, and leads 
to giant valley splitting in MoTe 2 . In order to gain 
further insight into the proximity-induced mag¬ 
netic interactions, a low-energy effective Hamilto¬ 
nian (LEH) is constructed to fit the band structure 
of MoTe 2 /EuO in the vicinity of K±, showing that 
Zeeman and Rashba fields are the key interactions 
induced by the proximity effect. The LEH is com¬ 
posed of four parts, H = EIq + Efsoc + 77ex + TLr. 
The four terms correspond respectively to the 
orbital interactions, SOC-induced spin-splitting, 
proximity-induced exchange and Rashba interac¬ 
tions, given as 


Ho = -b ayPy) + ycr^; (la) 

Hsoc = TS2{Xca+Xya-)- (lb) 

Hn = XniTSyOx - Sj:ay). (Id) 


Here, the Pauli matrices and Oa (a = 0, x, y, z) 
refer to real spin and orbital pseudo spin, respec¬ 
tively, and a± = ^(cq ± ct^). p is the electronic 
momentum and up is the Fermi velocity. The spin 
splitting of the conduction and valence bands due 
to intrinsic SOC is determined by parameters Ac 
and At,, respectively. The effective mass, m, cor¬ 
responds to the band gap in the absence of SOC, 
which results mainly from the crystal-field split¬ 
ting between dj ,2 and {dxy, d^^-y^} of Mo. Be and 
By are effective Zeeman helds experienced by the 
conduction and valence bands of MoTe 2 , arising 
from the coupling with a magnetic substrate. Note 
that the low-energy band structure of free-standing 
monolayer MX 2 , as shown schematically in Fig. 
1(c), can be described by Hq -b Etsoc (Eqs. (la) 
and (lb) together). [9] 

iLex represents the Zeeman field induced by the 
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substrate, and produces a band structure schemat¬ 
ically shown in Fig. 1(d). The valley degeneracy is 
clearly broken by the Zeeman field, with = 

—2Bc/v, which is independent of the valley index 
r. This, however, is inconsistent with the DFT re¬ 
sults (see Table I), for which is dependent on 

the valley index. To account for the r-dependence 
of and considering the surface electric field 

in the (lll)-direction, we include the Rashba field 
embodied in Hu- [30] The Rashba effect further 
hybridizes the valence and conduction bands, and 
mixes the spin components, with which the spin is 
no longer a good quantum number. However, we 
can still designate the low-energy bands of MoTe 2 
monolayer with up/down spin, since the out-of- 
plane spin components dominates in the valley re¬ 
gion (Fig. S2 of SI). With the Rashba term is 

added, become valley-dependent, in agree¬ 

ment with first principles band structure. Based 
on this model (Eqs. (la)-(ld) together), we can 
deduce that A^|,([ — = 2{By — Be). From the 

values in Table I, the effective Zeeman fields for the 
conduction and valence bands are seen to differ by 
as much as 30-40 meV, justifying the use of two ef¬ 
fective Zeeman fields, Bf./v, which reflect different 
effective Lande g-factors for the Bloch states. 

Matching the model with the first principles 
band structure of MoTe 2 /EuO leads to a semi- 
quantitative clarification of the role of proximity- 
induced magnetic interactions. First, we obtain 
the parameters, vp = 3.75 x 10^ m/s, Ac = —13 
meV, A„ = 107 meV and m = 861 meV, in the low- 
energy model of free-standing MoTe 2 monolayer 
(i.e. Hq + Hsoc), by fitting the corresponding DFT 
band structure. With EuO substrate is added, we 
keep Xc/v unchanged and all other parameters in 
Eq. 1 are refined by a least-square fit to MoTe 2 - 
projected bands near K± in Eig. 2(c) (detailed 
in SI). The Fermi velocity is largely unaffected by 
the substrate, with vp = 3.95 x 10® m/s. The ef¬ 
fective mass is significantly renormalized to a value 
984 meV. The effective Zeeman fields are gigantic: 
Be = 206 meV, By = 170 meV, which translates 
to a magnetic field over 2937 tesla. For a compar¬ 
ison, the magnetic field produced by the magnetic 
moments in Eu ions 4.4 A above the EuO(lll) sur¬ 
face is only ~0.007 tesla, which is computed from 
the direct summation of the magnetic field pro¬ 
duced by magnetic dipoles of Eu ions (7 /r^) in a 
semi-infinite EuO crystal. The Rashba parameter 
is determined to be Ar = 72 meV. The relevant en¬ 
ergy scales of the LEH model with as-determined 
parameters above are summarized in Table I, and 
it is seen that they are in decent agreement with 
the DFT results. It can be concluded that the 
proximity-induced spin interaction is dominated 
by gigantic Zeeman effect, with significant Rashba 
interaction. 

Discussions . — The valley splitting by 
proximity-induced magnetic interaction is at¬ 
tractive, as it creates giant difference in various 
energy scales between the valleys. In this section 


we discuss how the giant valley splitting will 
facilitate the access and manipulation of valley 
and spin degrees of freedom in a MoTe 2 /EuO 
heterostructure. In order to access the valleys of 
MoTe 2 , we will need to tune the chemical potential 
into the MoTe 2 gap (between ~ —0.9 and —0.4 
eV), which may be achieved with electrical gating 
or chemical doping. In the ensuing discussions, 
we will assume that we begin with a system of 
which the chemical potential is already tuned 
into the MoTe 2 gap. In what follows, we will 
discuss essential experimental consequences of the 
proposed valley degeneracy splitting, in terms of 
transport and optical properties, and tunability of 
valley splitting. 
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FIG. 3: Anomalous valley transport of hole-doped 
MoTe 2 /EuO heterostructure, (a) and (b) show the 
schematic hole-doped MoTe 2 -projected band structure 
of the heterostructure when the substrate EuO is mag¬ 
netized upward and downward, respectively (c./. Eigs. 
2 (c) and (d)). (c) Computed non-Abelian Berry cur¬ 
vature [26], D (units A ), of the valence bands occupied 
up to the MoTe 2 gap of the heterostructure when EuO 
substrate is magnetized upward (band structure in Fig. 
2(c)). The Berry curvature is momentum-resolved and 
shown in the entire Brillouin zone, (d) Schematic de¬ 
piction of anomalous Hall effects. The positive sign in 
circles denotes the hole. Holes with red up and blue 
down arrows denote up-spin one from K- valley in (a) 
and down-spin one from K+ valley in (b), respectively. 
E is the applied in-plane electric field. 

Owing to valley splitting, we can selectively 
dope only one valley with carriers at equilibrium. 
Based on the band structure in Fig. 2(c), it is 
seen that hole doping will give simple access to the 
valley states of MoTe 2 monolayer, whereas elec¬ 
tron doping will be interfered by the substrate’s 
bands. We will therefore focus on the case with 
hole doping. If we dope the heterostructure with 
holes at K+ valley when the substrate is magne¬ 
tized upward (Fig. 3(a)), the spin holes may pro¬ 
duce a transversal current under a longitudinal in¬ 
plane electric field (Fig. 3(d)). [31] The anomalous 








5 


transport arises from the intrinsic Hall conductiv¬ 
ity without magnetic held 

g2 r ^^2 

^xy ~ ^ J (^ 27 r)^ fl{k)n{k), (2) 

where il{k) is the Berry curvature of Bloch elec¬ 
tron [31], and n{k) the occupation of fc-space (from 
doping or otherwise), e the elementary charge and 
h the reduced Planck’s constant. This anomalous 
Hall conductivity is an intrinsic property of the 
Bloch band, and is the key to the detection of val¬ 
leys by electric measurement [1, 4], It therefore 
is desirable to have sizable B(fc) especially in the 
vicinity of valleys of of MoTe 2 monolayer. 

As shown in Fig. 3(c), the calculated Berry 
curvature is sharply peaked in the valley region, 
with opposite signs for K±. Clearly, the valley 
identity remains intact apart from the giant valley 
splitting and shall lead a transverse current un¬ 
der an applied electric held. The hux of the spin 
holes carries three observable physical quantities, 
namely, charge, spin moment and valley-dependent 
orbital magnetic moment, which can be referred to 
respectively as anomalous charge, spin and valley 
Hall effects. Due to the dominance of the Zeeman 
effect as discussed above, it is expected that when 
the magnetization of the substrate is hipped, the 
magnetically induced valley polarization will cor¬ 
respondingly be hipped. This conhrmed by the 
DFT calculation, as shown by the band structure 
in Fig. 2(d). Experimentally, this corresponds to 
switching the magnetization of EuO substrate with 
a relatively small magnetic held, as EuO has very 
weak magnetic anisotropy [32]. With the reversed 
magnetization, the valley polarization as well as 
the associated anomalous charge/spin/valley Hall 
effects are reversed (Figs. 3(b) and (d)). 

The valley identity is also associated with 
valley-contrasting circular dichroism in optical ab¬ 
sorption in addition to the contrasting anomalous 
transport, which is a very important property of 
MX 2 monolayer. The circular dichroism is readily 
characterized by the optical oscillator strength un¬ 
der the circularly polarized optical held, f±{k) = 
|7^±(fc)p/me(A/,fc - Ei^k), where V± = VxE iVy 
with Va = (■*/’/ba IV’i) (« = x,y) being the inter¬ 
band matrix element, mg is the electron mass, Ei 
and Ef are the energies of initial and hnal Bloch 
states, and ipi, ipj are the corresponding wavefunc- 
tions. As shown in Fig. 4(a) and (b), despite the 
giant splitting of valley degeneracy the optical ab¬ 
sorption of the MoTe 2 /EuO still preserves perfect 
circular dichroism in the vicinity of K± valleys. 
The optical oscillator strength for left-polarized 
light is sharply peaked near , but is vanishingly 
small near K-, and vice versa for right-polarized 
light, which will allow easy chiral optical pump¬ 
ing induced valley polarization in EuO-supported 
MoTe 2 monolayer, similar to free-standing MoTe 2 
monolayer. 

The proximity-induced Zeeman and Rashba ef¬ 
fects together lifts the optical frequency degener¬ 


acy, as discussed above, which adds an interesting 
variable to the dynamical optical valley polariza¬ 
tion. For free-standing MoTe 2 monolayer with de¬ 
generate valleys, non-polarized and linearly polar¬ 
ized light excite carriers equally from both valleys, 
and cannot lead to valley polarization. In contrast, 
dynamical valley polarization can be achieved in 
EuO-supported MoTe 2 monolayer by selection rule 
based on energy conservation, with a photon en¬ 
ergy huj satisfying the relation A[]pj < Hlo < A“pt 
for the EuO substrate magnetized upward, regard¬ 
less of light polarization. Due to the giant splitting 
of valley degeneracy, A^pj — A[Jpt = 44 meV, the 
energy selection works over a wide spectral range. 
When light within the energy range above illumi¬ 
nates the sample, only electron-hole pairs from the 
Ar+ valley are generated (Fig. 4 (c)). The photoex- 
cited electrons and holes may recombine, leading 
to only left-polarized luminescence. Under an in¬ 
plane electric field, the electrons and holes acquire 
opposite transversal velocities and move to oppo¬ 
site boundaries of the sample (Fig. 4 (d)), due to 
the opposite signs of the Berry curvatures in the 
conduction and valence bands. This also leads to a 
net charge/spin/valley Hall current. With the re¬ 
versed magnetization of the substrate, the carriers 
will be excited at K_ valley and the corresponding 
charge/spin/valley Hall effects are reversed. 



FIG. 4: Dynamical valley polarization in MoTe 2 /EuO 
heterostructure, (a) and (b) show the fc-resolved oscil¬ 
lator strength between the up-spin valence and conduc¬ 
tion bands, f{k), throughout the Brillouin zone for left- 
and right-polarized lights, respectively. The oscillator 
strength between the down-spin valence and conduc¬ 
tion bands is similar and not shown, (c) Schematic de¬ 
piction of optical pumping in MoTe 2 /EuO heterostruc¬ 
ture with EuO magnetized toward MoTe 2 monolayer. 
The lowest-energy interband transition of the two val¬ 
leys with energies of and A“pt are indicated, (d) 
Optically induced Hall effects. The positive and neg¬ 
ative signs in circles denote the photoexcited hole and 
electron, respectively. Up and down arrow denotes up 
spin and down spin respectively. E is the applied in¬ 
plane electric field. 
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The last and very important remark on the pro¬ 
posed valley splitting is concerning its tunability. 
As noted before, the magnetic anisotropy of EuO 
is weak [32] , and can be easily rotated with a rel¬ 
atively small external magnetic field. As this is a 
readily accessible experimental knob, it is useful 
to discuss how the valley splitting depends on the 
direction of substrate magnetization. As discussed 
above, the proximity-induced interaction is domi¬ 
nated by the Zeeman field, and we therefore ignore 
the change to the Rashba field as the magnetiza¬ 
tion rotates in a first approximation. A straightfor¬ 
ward generalization of Eq. (Ic) can be introduced 
to describe the field-tuning of substrate magneti¬ 
zation. 


— s * T ^ (3) 

where h is a unit vector denoting the direction of 
the proximity-induced Zeeman fields, which enters 
as Bc/v = nBc/v One notes that B^/v are the 
induced effective Zeeman field, rather than the ap¬ 
plied magnetic field for rotating the magnetization 
of the substrate. The latter field is supposed to be 
small so as not to interfere with subsequent trans¬ 
port measurements. 
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FIG. 5: Highly tunable valley splitting in MoTe 2 /EuO 
heterostructure, (a) Band structure in the vicinity of 
the valleys (momentum range 7r/5a around K±, where 
a is the lattice constant of MoTe 2 ) from the LEH 
model, for 9 — 7r/2,7r/4, 0, —7r/4, —7r/2. The bands 
are colored red and blue when the spins are antiparal¬ 
lel and parallel with the substrate magnetization (see 
insets), respectively, (b) Valley splittings as functions 
of the substrate magnetization direction. 

Taking fi = (cos 0,0, sind) and assuming all 
parameters of the model to remain unchanged, 
the corresponding evolution of MoTe 2 valley band 
structure is shown in Fig. 5(a) to demonstrate the 
effects of field tuning of substrate magnetization. 


It is evident that as the magnetization of the sub¬ 
strate is rotated, the valley splitting can be con¬ 
tinuously tuned. Compared to conduction bands, 
valence bands are more sensitive to the change of 
the magnetization direction. In Fig. 5(b), the val¬ 
ley splittings (see caption) for conduction and va¬ 
lence bands are plotted as functions of 6, as the 
magnetization rotates from upward, to horizontal, 
and finally to downward directions. As the mag¬ 
netic field turns from perpendicular to horizontal 
direction, the change in valley splitting for valence 
bands is ~ 46 — 67 meV, whereas the counterpart 
for conduction bands is ^ 3 — 19 meV. The large 
change of valence valley splitting is clearly advan¬ 
tageous, as the valence bands in the valleys are 
unobstructed by substrate bands or the bands in 
other regions of the A:-space (see Figs. 2 (c) and 

(d)). 

Conclusion . — In summary, we propose a gen¬ 
eral strategy to lift valley degeneracy in MX 2 com¬ 
pounds. As exemplified by computational model¬ 
ing of MoTe 2 /EuO heterostructure, this approach 
has several advantages. First, the valley split¬ 
ting is giant as demonstrated computationally for 
the MoTe 2 /EuO heterostructure, and much larger 
than the band shifts of ~ 0.1 meV/tesla produced 
by an external magnetic field. [15-18] Second, with 
the giant valley splitting, the valley and spin de¬ 
gree of freedom in MX 2 can be manipulated stat¬ 
ically and dynamically. Third, the induced valley 
splitting is highly tunable, by controlling the mag¬ 
netization of the EuO substrate. This is a clearly 
advantage compared to inducing valley splitting 
by magnetic doping [33]. Although we have dis¬ 
cussed only the valley splitting for the most sta¬ 
ble configuration of the heterostructure, the valley 
splitting is in fact also prominent for metastable 
configurations we have examined, as detailed in 
SI. The giant and tunable valley splitting adds 
a readily accessible dimension to the valley-spin 
physics with rich and interesting experimental con¬ 
sequences, and offers a practical avenue for explor¬ 
ing device paradigms based on the intrinsic degrees 
of freedom of electrons. On a final note, it should 
be emphasized that the effective Zeeman field ex¬ 
perienced by MoTe 2 is substantially enhanced by 
the interfacial interaction. The bare magnetic field 
on the EuO(Ill) surface is very small (~ 0.007 
tesla within the MoTe 2 ), which is shall not lead 
to significant quantum Hall effect with reasonable 
doping level at typical experimental temperatures. 
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